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Abstract In this paper, we study the weighted difference substitutions from geometrical views. First, we 
give the geometric meanings of the weighted difference substitutions, and introduce the concept of conver- 
gence of the sequence of substitution sets. Then it is proven that the sequence of the successive weighted 
difference substitution sets is convergent. Based on the convergence of the sequence of the successive 
weighted difference sets, a new, simpler method to prove that if the form F is positive definite on T„, 
then the sequence of sets {SDS ( - m ^(i r ')}^ =1 is positively terminating is presented, which is different from 
the one given in [TJ]. That is, we can decide the nonnegativity of a positive definite form by successively 
running the weighted difference substitutions finite times. Finally, an algorithm for deciding an indefinite 
form with a counter-example is obtained, and some examples are listed by using the obtained algorithm. 
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1 Introduction 

Theories and methods of nonnegative polynomials have been widely used in robust control, 
non-linear control and non-convex optimization [l|, 0, etc. Some famous research works on 
nonnegativity decision of polynomials without cell-decomposition were given by Polya's Theorem 
0,(1 and papers @,0j. 

Few years ago, Yang H & OH 

introduced a heuristic method for nonnegativity decision of 
polynomials, which is now called Successive Difference Substitution (SDS). It has been applied 
to prove a great many polynomial inequalities with more variables and higher degrees. Yang 
recommended further studies on SDS and put forward some open problems. 

A valuable progress on the topic was made by Yao [ll[ . He investigated a weighted difference 
substitution 

Xi = h + — -\ 1 , 

2 n 

t% , tn 

X2 = — H 1 , 

< 2 n (1) 



%n — ? 

k n 

instead of the original difference substitution 
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x\ = ti + 1% H V t 

x 2 = t 2 + h t n , 



Xn tn j 

and proved that, for a form (namely, a homogeneous polynomial) which is positive definite on R™ , 
the corresponding sequence of SDS sets is positively terminating, where R™ = {(x±, x%, x n )\xi > 
0, i = 1, 2, n}. That is, we can decide the nonnegativity of a positive definite form by successively 
running SDS finite times. 

This paper is organized as follows. Section 2 introduces some preliminary notions of the 
weighted difference substitutions. Section 3 provides a new perspective to study the weighted 
difference substitutions, gives the geometric meanings of them, and proves that the sequence of 
the successive weighted difference substitution sets is convergent. A new, simpler method to 
prove that if the form F is positive definite on R™, then the sequence of sets {SDS^ m ^(F)}^ =1 is 
positively terminating is presented in Section 4, and an algorithm for deciding an indefinite form 
with a counter-example in Section 5. By using the obtained algorithm, several examples are listed 
in Section 6. 



2 Preliminary notions 



We first introduce some notations and definitions according to (with some differences). 
Consider W n G R nx ™, where 



W n = 



1 







? 



o h 



(3) 



Let [k% k% ■ ■ ■ k n ] be a permutation of 1,2, ••• ,n. P[kik 2 ---k n ] = [ a (*jj)] is an n x n matrix 
for which a(l,ki) = l,a(2,&2) = , o(n, fe n ) = 1, and in all other positions (Permutation 

matrix) . 

Definition 2.1. n x n square matrix B^ 1 k 2 ---k n ] is defined as follows: 

B[k 1 k 2 —k n ] = P[k 1 k 2 ...k n ]W n . 

And the set that consists of all -B[fc 1 fc 2 ...fc„] is denoted by PW n (there are n! elements in PW n ) and 
called the weighted difference substitution matrix set. Accordingly, the set of linear transforma- 
tions 

{X Tr = B [a] T Tl \B [a] ePW n }, 

is called the weighted difference substitution set, which consists of nl substitutions, where X, T G 
R" , and X Tl , T Tr are respectively the transposes of A, T. 

Definition 2.2. The set of linear transformations 

{A Tl = B [ai] B [a2] ■ ■■B [am] T^\B [ai] G PW n }, 

is called the m-times successive weighted difference substitution set, which consists of (n\) m sub- 
stitutions. 

Definition 2.3. Given the form F G R[xi,X2, ■ ■ ■ ,x n ], when [ai], [02], • • ■ , [a m ] traverse all the 
permutations of 1, 2, • • • , n respectively, we define the set 



n! n! n! 



SDS^(F)= |J - ..(J \jF(B [ai] B [a2 yB [am] X T <), 

[ot m ] l a 2] [oil] 

which is called the m-times successive weighted difference substitution set of the form F. 



Definition 2.4. We define the sequence of sets {SDS (m) (F)}^ =1 as follows 

{SDS(F) (m) }~ =1 = SDS(F),SDS (2) (F),--- . 

It's time to define the termination of {SDS^^F)}^—!, which is directly related to the positive 
semi-definite property of the form F. 

Let a — (oti, &2, • • • , ot n ) G N™, and let \a\ = ai + a 2 + • • ■ a„. Then we write a form F with 
degree d as 

r > C a Xi X2 ' ' ' X n . 

\a\=d 

Definition 2.5. The form F is called trivially positive if the coefficients c Q of every term 
x^x^ 2 ■ ■ ■ x" n in F are nonnegative. If F(l, 1, • ■ ■ , 1) < 0, then F is called trivially negative. 

Definition 2.6. If the form F(X) > for all X G R™, then F is called positive semi-defmite ; 
If F(X) > for all Xty 0) G R™, F is called positive definite on KTf.; If there are X and Y G E™ 
such that F(X) > and < 0, F is called indefinite on K7J.. 

Lemma 2.1. Given the form F on R™, if the form F is trivially positive, then F is positive 
semi-definite on R™ ; If the form F is trivially negative, then F isn't positive semi-definite on R" . 

Definition 2.7. Given a form F on R™ , if there is a positive integer k such that every element of 
the set SDS (fc) (F) is trivially positive, the sequence of sets {SDS (m) (F)}~ =1 is called positively 
terminating; If there is a positive integer k and a form G such that G G SDS^(F) and G is 
trivially negative, the sequence of sets {SDS'"^ (F)}^ =1 is called negatively terminating; The 
sequence of sets {SDS^ m '(F)}^ =1 is neither positively terminating nor negatively terminating, 
then it is called not terminating. 

By Definition 2.7, it's easy to get the following lemma. 

Lemma 2.2. Given the form F on R™, if the sequence of sets {SDS (m) (F)}^ =1 is positively 
terminating, then F is positive semi-definite on R"; If the sequence of sets {SDS^ m ^(F)}^ =1 is 
negatively terminating, then F isn't positive semi-defmite on R™ . 

Next, we will give the definition of the normalized substitution. 

n 

Definition 2.8. Let V = [uy] be an n x n matrix. If v ij = 1, J = 1,2, ••■n, V is called a 

i=l 

normalized matrix. And the corresponding substitution 



X Tr = y T Tr 



is called a normalized substitution. 



Lemma 2.3. Let U = [uy] = V1V2 • ■ ■ Vk- If Vj., V2, • ■ ■ ,14 are normalized matrices, then U is a 
normalized matrix , that is, 

n 

^2 u i] = l(j = 1,2, • • -n). 

i=l 

n n 

And let (a?i, X2, • • • , £n) Tr = • ■ • , tn)' 11 , then Xi = 1 iff ti = 1. 



The proof of Lemma l2.3l is very straightforward and is omitted. 
The (n — l)-dimensional simplex is defined as follows 



T„ = {(xi,X2,--- ,x n )\y]xj = l,(xi,x 2 ,--- ,x n ) el"}. 

i=l 



Definition 2.9. If the form F(X) > for all X G T n , then F is called positive semi-definite 
on T„; If F(X) > for all X G T n , then F is called positive definite on T„; If there are X and 
Y ET n such that F(X) > and F(Y) < 0, then F is called indefinite on T n . 

Obviously, we have the following conclusion. 

Lemma 2.4. The form F is positive semi-definite (positive definite, indefinite) on T„ iff F is 
positive semi-definite (positive definite, indefinite) on K™ . 

According to Lemma 12.41 for brevity, we suppose that the form F is defined on T n in the 
remainder of this paper. 

3 Convergence of the sequence of successive weighted difference 
substitution sets. 

In this section, we'll consider the weighted difference substitutions from geometrical views. 
Let X = (xi, X2, ■ ■ ■ ,x n ) G T„. Consider the weighted difference substitution 

X T * = W n T T \ (4) 

where W n is denoted by ©. 

By ©, if X = Ax = (1,0,. ..,0), then T = (1,0,...,0); If X = A % = (|,|,...,0), then 
T = (0,1,...,0); If X = A n = (i i,...,!), then T = (0,0,...,1). Moreover, A k is the 
barycenter of the (k-l)-dimentional proper face of T„ which contains the points A\, A^, ■ ■• , Ak for 
k = 1, 2, • • • , n . According to Lemma l2.31 for all (a?i, x^, ■ ■ ■ , x n ) € T„, we have t\ + ta + ■ • • + £«, = 
1. Therefore, A1A2 ■ ■ ■ A n is a subsimplex of the first barycentric subdivision of T n , satisfying 
W n = [Aj r , A2 1 , ■ ■ ■ jAj^]. And it is indicated that the weighted difference substitution 0$ and 
the subsimplex A1A2 ■ • ■ A n correspond to each other. 

Analogously, the other nl — 1 weighted difference substitutions correspond to the other nl — 1 
subsimplexes of the first barycentric subdivision of T n . 

Hence, from geometrical views, the weighted difference substitution set corresponds to the first 
barycentric subdivision of T„, that is, there is a one-to-one correspondence between the weighted 
difference substitutions and the subsimplexes of the first barycentric subdivision of T„ . 




Figure 1: Barycentric subdivision 



For instance, when n — 3, the following six weighted difference substitution matrices correspond 
to subsimplexes labeled as 1-6 in Fig. 1, respectively, 
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1 I 
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By Lemma 12.31 and Definition 12.21 we have that the m-times successive weighted difference 
substitution set corresponds to the m-th barycentric subdivision of T n . 

Next, we'll introduce the concept of convergence of the sequence of substitution sets. 



Definition 3.1. Let a be a subsimplex of T„, the maximum distance between vertexs of a is 
called the diameter of a. 



Definition 3.2. Let Kq — T„, and Ki + \ be the subdivision of if, for i = 1, 2, • • • . If for all 
e > 0, there exists TV G N such that all the diameters of the subsimplexes of Km are less than e, 
the subdivision sequence {Ki}°l 1 is called convergent. 

Definition 3.3. Suppose that the subdivision scheme through which is subdivided into Ki 

corresponds to the substitution set Li for i = 1,2,---. If the sequence {-fQli^i is convergent, the 
sequence of substitution sets {L i }°°^ 1 is called convergent. If Li = L, i = 1, 2, • • • , briefly, we say 
the sequence of the successive L-substitution sets is convergent. 

It's time to consider the convergence of the sequence of the successive weighted difference 
substitution sets. 

Lemma 3.1. 12, HI Let K be a complex. If Kn is the fc-th barycentric subdivision of K, then 
for all e > 0, there exists N € N such that all the diameters of the subsimplexes of Km are less 
than s. 

By Lemma 13.11 , we have the following theorem, which plays important roles in the proofs of 
Theorem 14 . 1 1 and Theorem 15. II in the next sections. 

Theorem 3.1. The barycentric subdivision sequence {K m }^ =1 of T n and the corresponding 
sequence of the successive weighted difference substitution sets {SDS*-™ 1 ' {F)}™ =1 are convergent. 



4 Nonnegativity decision of forms 

Given a form F on T n . We know that if the sequence of sets {SDS^™- 1 (F)}™ =1 is negatively 
terminating, then we can conclude that F is positive semi-definite on T n . Thus there is a natural 
question, that is, which kind of forms can be solved by the method? Yao [ll[ proved that, for a 
positive definite form F on T n , the sequence of sets {SDS^(-F) is positively terminating. 

In this section, we present a new, simpler method to prove the conclusion, which is based on 
the convergence of the sequence of the successive weighted difference substitution sets. 

Theorem 4.1. Let the form F be positive definite on T n , then the sequence of sets {SDS^ (F)}™_ 
is positively terminating. 

Proof. We only give the proof for the ternary form with degree d, and the multivariate form can 
be gotten by induction. 

Suppose that F(x%, X2,Xs) = aijkx\x\x\. An arbitrary m-times successive weighted 

i-\-j-\-h=d 

difference substitution can be written as 

x\ = k\t\ + (ki + a 1 )t 2 + (ki + /3i)i 3 , 

x 2 = kah + (fca + a 2 )t2 + (k 2 + fo)h, (5) 

2>3 = Ml + (&3 + «3)*2 + (&3 + ^3)*3, 
3 3 3 

where fc, = 1, ^ a, = and /3j = 0. 

i— 1 i— 1 i— 1 

Let t = t\ + t 2 + ts, then §5§ becomes 



X\ = kit + ait 2 + ftit 3 , 
x 2 = k 2 t + a 2 t 2 + p 2 t 3 , 
X3 = fat + a 3 t 2 + P 3 t 3 . 



(G) 



Thus 

$(*i,*2, t 3 ) = F(kit + aite + ^1*3, M + a 2 t 2 + f3 2 t 3 , k 3 t + a 3 t 2 + (3 3 t 3 ) 

= £ ^3k{kit + ai t 2 + (3 1 t 3 ) l (k 2 t + a 2 t 2 + (3 2 t 3 y(k 3 t + a 3 t 2 + p 3 t 3 ) k 

i-\-j-\-k—d 

p+q+r=j,p=ij p+q+r=k,p=ik 



= £ a ljk klk 3 2 k 3 t d + £ <p ljk {ki,k 2 ,k 3 ,ai,a 2 ,a 3 ,l3i,(3 2 ,l3 3 )t\t 3 2 t 

i-\-j-\-k—d i-\-j-\-k=d 

= F{k u k 2 ,k 3 )t d + £ <^ijfc(fei,fci)A;3,Q!i,a2,a3,/3i,/92,/93)*i*2*3 
= £ f-y^^(fci,fe, fc 3 ) + ^ijfe(fci,/c 2 , fc3,ai,a2,a3,/?i,^2,/33) ) *1^*3 

i+j+k=d 

= £ Aijkt\t{t 3 . 

i+j+k—d 



( 7 ) 



where 



-^ij-A = Tpj^^(fci,fc2,fc 3 ) + 0. iife (fci,fc2,/c3,ai,a2,a3,/3i,/?2,/33)- (8) 

Obviously, 

lim <fcj*(fei,fe,fc3,ai,a2,a3}/3i)/32,/33) = 0. fQ\ 

(ai,aa,a3.^i,/32 ) /9 3 )->(0,0,0,0,0 1 0) V ' 

And since -F(a^i, x 2 , x 3 ) is positive definite on T„, there exists e > such that 

k 2 , k 3 ) > £ > 0. (10) 

On the one hand, the vertexs of the subsimplex which corresponds to the successive weighted 
difference substitution ([5]) are respectively 

(ki,k 2 , k 3 ), (fci + ai,fe 2 + a 2 ,k 3 + "3), (hi + Pi,k 2 + /3 2 , k 3 + (i 3 ). 

By Theorem (|3.1jl . ai, a 2 , a 3 , /3i, /3 3 can be sufficiently small when m is sufficiently large. 

On the other hand, for F is continuous on T n and by (j5j)- (fTU|) , we have Aijk > when 
ai, a 2 , a 3 , f3i, p 2 , j3 3 are sufficiently small. 

Putting together the above two aspects, we have that there exists a sufficiently large integer m 
such that F becomes trivially positive by (|5|). For the successive weighted difference substitution 
([5]) is arbitrary, the sequence of sets {SDS^™- 1 (F)}™ =1 is positively terminating. □ 

According to the proof of Theorem 14. 11 we obtain the following conclusion. 

Corollary 4.1. Let the form F be positive definite on T„, then by an arbitrary m-times succes- 
sive weighted difference substitution, when m is sufficiently larger, F can become a nonlacunary 
trivially positive form. 

Theorem 14. II is somewhat analogous to Polya's Theorem. However, many examples show that 
Polya's Theorem seems almost useless to positive semi-definite forms except for few cases, while 
SDS is demonstrated very helpful to positive semi-definite ones as well. 



5 Decision of indefinite forms 



Many problems, such as the inequality disproving, are always transformed into decision of 
indefinite forms. 

Given a form F on T n . Suppose that there exists Xq G T„, such that F(Xo) > 0. It is 
well-known to us that if the sequence of sets {SDS < - m - ) (i' 1 )}^ =1 is negatively terminating, then F 
is indefinite on T„. Then it follows a question naturally: for an indefinite form F on T„, is the 
{SDS' m - > (i r )}m = i negatively terminating? The following theorem answers the question. 

Theorem 5.1. Let the form F be indefinite on T„, then the sequence of sets {SDS^ 1 ^)}^! 
is negatively terminating. 

Proof. Since the form F is indefinite on T„, there exists X$ G T„ such that F(Xo) < . And F 
is continuous on T„, so there exists a neighborhood IJ(Xq) C T n of Xo (If Xq is on the boundary 
of T„, then we take V(X a )f]T n ) such that F(X) < for all X G U(X ). For the barycentric 
subdivision sequence of T„ is convergent, then there exists a subsimplex a of the fc-th barycentric 
subdivision of T„ , which corresponds to the fc-times successive weighted difference substitution 

* Tr = B [tl] B [l2] ■ ■ ■ B [ik] T^, • • , B [lk] G PW n , 

satisfying a C U(Ao), where A: is a sufficiently larger integer, and PW n is the weighted difference 
substitution matrix set. Thus, —F(X) is positive definite on a. By Theorem 14.11 the sequence 
of sets {SDS*-™- 1 (— F)}m = i is positively terminating, so there exists an /-times successive weighted 
difference substitution 

xTl = B [ii\ B [h\ ■ ■ ■ B l:ii} TTl > B IJi] > B [h\ > ' ' * > U Ji. e PW n, 

satisfying that F{B^B^ ■ ■ ■ Bu^BsaaB^i ■ ■ ■ S[ Ji ]T Tr ) is trivially negative. Therefore, the se- 
quence of sets {SDS^^F)}^]^ is negatively terminating. □ 

By the proving process of Theorem 15.11 we obtain the following algorithm, which is used to 
decide the nonnegativity of a form or to decide an indefinite form with a counter-example. 
Algorithm (YYS) 

Input: the form F £ <Q[x\,X2, ■ ■ ■ , x n ], where F is positive definite or indefinite on R™ . 
Output: "The form F is positive semi-definite", or "Xo, F(X ) < 0". 
stepl: Let F = {F}. 

step2: Compute (J SDS(F), 

Few 

Let 

F = (J SBS{F) - {trivially positive forms in [j SDS(T)} = {F {1] , F [2] , • • • , F [k] }, 
fsf FeF 

where 

F[i\ = F(B {i] X^), B {i] g PW n . 

step3: Let L= [[1],[2],--- , [k]]. 

step31: If F is null, then output "the form F is positive semi-definite", and terminate. 
step32: If there is a trivially negative form Fu] G F, then output 

"X = B m (-,-,...,-) Tl , F(Xo) <0", 
1 ' n n n 

and terminate, where 



B L[i] - B [L[m] B [L[m] ■ ■ ■ B [m[m]] , 



L[i] is the z-th component of L, L[i]\j] is the j-th component of L[i], and m is the 
the total number of the components of L[i]. 
step33: Else, Compute [j SDS(F). Let 

Few 

F = |J SDS(F) - {trivially positive forms in (J SDS(F)} 
= {^[o P (L[i]),i] , ■ • • , F [op{L[1])M] , F [op(L[2 ])4] , • • • , F [op{L[2])M] , 

■■■ ,-F[op(£[fc]),l]> • • ■ , F [op(L[k]),l k ]}, 

where 

*[op(L[i]),j] = ^(-B[op(L[i]),j]^ Tr ), 

and op(L[i]) extracts operands from op(L[z]). And let 

L = [[op(L[l]), 1], • ■ ■ , [op(L[l]), ii], [op(L[2]), 1], • ■ ■ , [op(L[2]), h], 
■■■ ,[op(L[fe]),l],--- ,[op(L[k]),l k }}, 

then go to step3. 

By Algorithm YYS, we design a Maple program called YYS, see Appendix. To the program 
YYS, there are some positive semi-definite forms making the program do not terminate, that is, 
we cann't decide these positive semi-definite forms by the method. 

6 Examples 

In this section, we demonstrate the program YYS with some examples. 
Example 1. Show that the following form is positive semi-definite on R^_, 

F(x, y, z) = x(x - yf - y(-z - yf - z(x - zf. (11) 

Utilize the program YYS and execute order YYS(F,[x,y,z]). The procedure need successively 
run the weighted difference substitutions 3 times, then outputs: "The form F is positive semi- 
definite" . 

Example 2. Show that the following form is indefinite on R+. 

F(x, y, z) = 7x 3 - I2x 2 y - 12x 2 z + 6xy 2 + \2xyz + 6xz 2 - ^y 3 - 3y 2 z - 3yz 2 - ^z 3 . (12) 

Executing order YYS(F, [x,y,z]), we have a counter-example by successively running the 
weighted difference substitutions 2 times: 

A = (^,i^f, F(A> )<0. 
v 108' 108' 54 7 ' v °' 

Obviously, F(l, 0, 0) = 7 > 0, so the form F is indefinite on R^_. 

Example 3. Let x>0,y>0,z>0, and x + y + z ^ 0. Try to decide whether the following 
inequality holds for all p e N. 

3 y + z z + x x + y' y 3 ' ~ y ' 

Take off denominators of the left polynomial, and denote the new polynomials by F\,F2,- ■ ■ ,Fq 
for p = 1, 2, ■ • • , 6, respectively. Execute order YYS(F p , [x,y,z]), p = 1, 2, ...6. For p = 1, 2, • • • ,5, 
the procedure only need run the weighted difference substitutions 1 time, then outputs: "The 
form F is positive semi-definite" . For p = 6, we have a counter-example by successively running 
the weighted difference substitutions 5 times: 

- 2159 3685 3661 T - 

x °- Wneei'IIeeT ' F e(Ao)<o. 



So the inequality cann't hold for p = 6. 
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Appendix. Maple Program YYS 

YYS : =proc (poly , var) 

local a,b,A,f,i,j,k,m,n,r,s,t,F,G,H,M,W, newvar , st , Var : 

uses combinat, LinearAlgebra: 

F: = [[poly, [0]]] : 

n:=nops(var) : 

Var :=convert( var, Vector) : 

W:=(n)->Matrix(n,n, (i,j)->'if ' (i<=j , 1/j ,0)) : 
b:=permute(n) : a:=W(n) : 

A: = [seq(<seq(a[b[i] [j]] , j=l . .n)> , i=l . .n!)] : 
for i to nops(A) do 
for j to n do 

newvar [i , j] :=op( j , convert (A [i] . Var, list) ) : 

od: 

od: 

r:=100: 

for s to r do 

m:=nops(F) : 

f : = [] : 

for k to m do 
G[k]: = []: 

for i from 1 to nops(A) do 



st :={seq(Var [j] =newvar [i , j] , j=l. .n)}: 

G[k] : = [op(G[k]), [expand ( subs (st,F[k] [1])) , [op(F[k] [2]),i]]] : 

od: 

od: 

F: = [seq(op(G[u]) ,u=l. .nops(F))] : 
for i to nops(F) do 

if max([coeffs(F[i] [1])]) < then 
print (F[i] [2] ) : 
M:=IdentityMatrix(n) : 
for j from 2 to nops (F [i] [2] ) do 
M:=M.A[F[i] [2] [j]] : 

od: 

print (convert (M. Vector [column] (n, 1/n) , list) ) : 
return ("The form is indefinite") : 
elif min( [coeff s(F[i] [1] )] ) > then 
f : = [op(f),i] : 

fi: 

od: 

if nops(f)>0 then 

F:=subs({seq(F[f [t]]=NULL,t=l . .nops (f ) )},F) : 

fi: 

if nops(F)=0 then 
print (s) : 

returnO'The form is positive semi-definite"): 

fi: 

od: 

end proc : 



